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Abstract 

A phenomenological 2D model, simulating the martensitic transformation, is built upon existing experimental observations 
that the size of the formed plates -in direct transformation- decreases as the temperature is lowered; then they transform back in 
reversed order. As such, if a reverse transformation is incomplete (’’arrested”), the subsequent direct one will show anomalously 
large number of big size plates-old plus newly formed- but consequentially a depletion of intermediate sizes, due to geometrical 
constraints, phenomenon that generates thermal memory. 
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1. Introduction 

The temperature memory effect, as manifested by shape 
memory alloys, remarkably shows that a thermodynamic sys¬ 
tem can remember not only the sense of a transformation, 
as in magnetic hysteresis, but also quasi-precise temperatures 
at which certain actions were performed in the past, such as 
the interruption of a previous phase transition. Naturally, the 
subject attracted a number of experimental efforts (see, e.g. 
GlEIElllllIllSlIIlIlllSlIIol towards a clear description, but sur¬ 
prisingly few theoretical scenarios have been proposed (see e.g. 
Gnaa ), an unanimously accepted phenomenological model 
being so far rather elusive. 

Succinctly, the temperature memory effect (also called ’’ther¬ 
mal arrest”) consists in the system remembering one -or 
several- previous incomplete phase transition(s), as described 
in the following four-steps process: at step I, a complete direct 
phase transition (parent to product) is performed, then at step 
II the reverse (product to parent) transformation is ’’arrested” - 
stopped before completion- at a temperature ’’Ta”. Next, at step 
III, another complete direct transformation is performed and fi¬ 
nally at step IV a complete reverse one. In this final reverse 
transformation, the calorimetric signal shows a dip at a temper¬ 
ature close to Ta, at which the transformation was previously 
arrested at step II. 

Efforts are on-going towards understanding of the thermal 
memory effect. Previous theoretical works suggest that the re¬ 
lease of stress from the untransformed martensite can increase 
its transformation back temperature, which results in a sepa¬ 
rated second peak EllTlEl. j .Rodriguez-Aseguinolaza et al. la 
propose a more quantitative model also relaying on the redistri¬ 
bution of stress amongst martensite fractions and the parent ma¬ 
trix after one or multiple incomplete reverse transformations. 

In this paper we propose a different mechanism possibly rel¬ 
evant for the thermal memory effect, namely that incomplete 


thermal circles influence the distribution of the product phase 
plates sizes, making the latter a ’’witness” of the thermal his¬ 
tory - which should be readable by a calorimetric scan. 

The idea is rather simple and is based on the existing micro¬ 
scopic evidence that the martensitic transformation takes place 
with formation of plates, which do not grow indefinitely, but 
rather stop at a certain maximum size; then, the phase transi¬ 
tion proceeds by formation of other plates. There is also evi¬ 
dence that this intrinsic maximum size decreases with tempera¬ 
ture (see, e.g. nil). Keeping in mind that the size of a formed 
plate possesses such an intrinsic limitation, it is a simple next 
step to assume that geometrical constraints may impose supple¬ 
mentary limitations: the new plates can be limited in growth by 
the existing puzzle of neighboring plates. As such, the distri¬ 
bution of plates sizes becomes history dependent and it will be 
shown that this property can generate thermal memory. 

The outline of the paper is as follows: in Section 2, the gen¬ 
eral assumptions of the model are described, then a numerical 
illustration is given in Sections 3 and 4; Section 5 discusses 
briefly the relevance and limitations of the calculations, and 
Section 5 concludes the paper. 

2. Model 

In this section we present the basic assumptions within a sim¬ 
plified phenomenological model for a solid state phase transi¬ 
tion inspired by certain properties of the martensitic transfor¬ 
mation. The basic features of the proposed model are: 

(i) Finite size plates. The product phase, once nucleated, grows 
in form of plates only up to some intrinsic maximum 
size where the growth stops. Then, the phase transition 
proceeds by formation of new plates (rather than further 
growth of existing ones). In the numerical simulations. 
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this intrinsic maximum size will also be considered to de¬ 
crease with temperature lIT^ . 

(ii) Key role of geometrical constraints. The plates may stop 
growing sooner if they encounter other -already existing- 
plates which get in the way of their growth. The geometri¬ 
cal constrains thus work in one sense, of reducing the size 
of certain plates, and in this way influencing the plates dis¬ 
tribution statistics. 

(iii) Reverse transformation in reversed sizes order. In the re¬ 
verse transformation, the plates return to the parent phase 
in the reverse order of their formation, i.e. the smaller 
plates which were the last to form are the first to transform 
back and the largest plates transform back last. 

The property (i) already indicates a structured product phase 
composed of plates of different sizes, which can in principle 
enclose information. 

Property (ii) imposes a dependence on the specific sample 
history meaning that the number and sizes distribution of pre¬ 
existing plates statistically influence the sizes of newly formed 
plates. This is the proposed key mechanism responsible for the 
thermal memory effect. 

The property (iii) describes the reverse transformation. Its 
justification relies in the higher surface-to-volume ratio for 
smaller plates which makes them the first to become thermody¬ 
namically unstable (a numerical exemplification will be given 
in Section 4). 

The above described phase transition is schematically illus¬ 
trated in Fig.l. On the upper row, a direct transformation is 
indicated, in which the larger ’’martensite” plates nucleate and 
grow first (the red squares in panel 2), then as temperature de¬ 
crease, intermediate plates form (yellow squares in panel 3), 
and Anally the smaller blue squares form at lowest temperatures 
(panel 4). 

Let us now assume that in the reverse transformation marten¬ 
site plates transform back in reversed order, starting the small¬ 
est blue squares and followed by the yellow ones (4-2 reverse 
sequences). Now if the reverse transformation is stopped (’’ar¬ 
rested”) before the bigger red plates can transform back, they 
will remain untransformed and a subsequent direct transfor¬ 
mation will start from this pre-existing plates distribution, as 
depicted in Fig.l, second row (panels 2-5-7), and one can no¬ 
tice the ’’anomalously” large number of big size plate (those 
remained untransformed and the newly formed ones). A first 
direct effect is seen in panel 6, corresponding to temperatures 
favoring the yellow intermediate plates, but which now cannot 
grow to their normal size because they cannot fit in the existing 
puzzle; instead the formed germs grow only to blue sizes. The 
final sizes distribution is considerably changed. 

The limitation of plates growth, their random positioning in 
the parent matrix, and the geometrical constrains in growth are 
non-equilibrium effects, reflecting the pronounced metastable 
character of the transformation. 

The intuitive scenario depicted in Fig.l shall be numerically 
tested in the following section. 
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Figure 1: Schematic representation of thermal memory encrypted in the dis¬ 
tribution of plates sizes. The main assumptions made are that the sizes of the 
formed plates decrease with temperature (T\ > T 2 { 5 ) > T^ie) > ^ 4 ( 7 )), and they 
transform back in reverse order of their formation. A first direct transforma¬ 
tion (1 ^ 4) is followed by an incomplete reverse one (4 ^ 2) and a second 
direct transformation (2 ^ 7). Sizes distribution for the panels 4 and 7 are 
indicated by the bar charts. In the (2 ^ 7) transformation, intermediate size 
plats (yellow) are missing, because they cannot fit in the puzzle defined by the 
’’anomalous” large number of bigger (red) plates, therefore they only grow to a 
’’blue” size. 

3. Numerical simulations assuming isothermal direct trans¬ 
formation. 

In our numerical simulation, the product phase will be as¬ 
sumed to nucleate in random positions in the parent matrix and 
to grow in the form of squares which gradually All the avail¬ 
able surface. If a square, during its growing process, encounters 
neighboring (already existing) squares, or the system borders, 
it will stop growing, remaining at that respective size. 

While the 2D model is justifled by the plates-like growth of 
martensite germs (also, the alloys can be obtained as thin Aims, 
etc.), choosing a square shape for the plates should be regarded 
as a simplifying assumption- with the goal to just offer an illus¬ 
tration of the phenomenology. 

It is a good point to mention that an interesting debate is still 
on-going whether the martensite transformation is isothermal or 
athermal cniiiEiiisiiniiia, the balance seemingly leaning 
towards the isothermal side, even if sometimes with an atypical 
long transformation time |[T9l|20l[^. We shall assume the sim¬ 
plest formula for a isothermal nucleation rate (denoted /a^m), 
governed by the activation critical free energy AQc(7') (corre¬ 
sponding to the formation of the critical cluster) 12211231 : 

JA^uiT) = JoS exp{ - -^^1 (1) 

For our 2D case, the free energy gain when a product phase 
germ is formed (which in our model is a square of side L) can be 
written: ADa^m = L^-AF-l-4L-(T.The condition ^AQa^m = 0 
yields the critical square side = -2 (t/AF, from which it re¬ 
sults: ADc(7") = -4(T^/AF, AF = e{T - Tq) being the free 
energy difference in the linear approximation. Tq is the equi¬ 
librium transition temperature; e and a are respectively the en¬ 
tropy difference (between parent and martensite phase) and the 
interface energy per unit length. The interface energy is always 
positive and may be a sum of elastic and inelastic contributions. 
The factor S in Eq. 1 is the untransformed remained surface (the 
initial total surface is considered ). /q is a proportionality 
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constant, its value choice for the numerical simulations being 
described below. 

The quantities introduced so far are typical for a first order 
phase transition, but now we need to impose the supplementary 
specific property that the formed germs have an intrinsic limi¬ 
tation in growth- moreover, their final size also decreases as the 
temperature is lowered. One possible way to introduce intrin¬ 
sic limitation in growth is to assume a proportionality between 
the size of the critical germ (L^) and the size of the final plate 
{Lf) :Lf = nLc [meaning Lf ~ 1/(7 - Tq)]. One can notice 
that the size Lc of the critical germ decreases with decreasing 
temperature below To, then it results from our proportionality 
hypothesis, that the final size of the plates decrease as well with 
temperature. The phenomenologically introduced limitation in 
the plates growth completes our model for the direct transfor¬ 
mation for which numerical results are presented in this section. 

In order to proceed with the real-time numerical simulation 
of the phase transition, one needs to describe also the heat 
exchange with the calorimeter, which is considered simply: 
dQjdt = - 7'), where 7 and 7' are, respectively the 

sample and thermostat temperatures and A} is the total sample 
surface. In the numerical simulations, the thermostat temper¬ 
ature will vary linearly, the typical case in DSC experiments: 
(say, 7' = at, with a < 0 for the direct transformation). If the 
phase transition is initiated, part of the transferred heat is used 
as latent heat (2 l), the rest warms (or cools) the sample with 
the amount per unit time: dTjdt = d(Q - QL)ldt • l/(cA^) , 
where c is the specific heat. Ql is calculated as follows: the 
germs formation ratio is given by Eq.l, then each germ has 
a size dependent growing time, let this time be denoted t{L)\ 
then dQildt = 2 l/[r(L)] • L^De, the summation being over all 
growing germs, De being the latent heat per surface unit. In 
the numerical calculations the border energy will also be con¬ 
sidered. 

A brief discussion is in order about the parameters used in the 
numerical simulation. Since our intention is a purely qualitative 
insight on the phenomenon, the parameters are adimensional, as 
used in the numerical simulation. First, the sample is modeled 
by a square of 120X120 (discrete points). Much like the finite 
differences method for solving differential equations, the dis¬ 
cretization brings an important advantage by reducing the com¬ 
puting time, with the draw-back that the continuous space is just 
approximately described. The product phase plates are squares 
with the side ranging from Luax = 12 points, to Lum = 3. 

The equilibrium transition temperature is generically cho¬ 
sen as 7o = 300. The parameters appearing in the expo¬ 
nential from Eq.l are chosen such as 4cr^KeksTo) ~ 700 , 
Jo = 0.1, which lead to a ’’martensite start” temperature of 
approx Tms = 7o-15. Technically, Tms can be calculated 
from the condition that the temporal integration of the nucle- 
ation rate reaches the value 1, corresponding to the formation 
of the first germ JA^MiT{t))dt = 1. Other parameters (de¬ 
fined in this section) were chosen such that (DEL^^^J/tiLMax) < 
PA^(To - Tms)^ meaning that the latent heat is in general less 
that the total heat exchanged with the reservoir so that the direct 
phase transition, for instance, is accompanied by a permanent 


cooling of the sample, making possible the formation of pro¬ 
gressively smaller plates. 

The first plates to form are also the largest, and in our model 
the maximum size will be Luax = 12, from which the propor¬ 
tionality constant assumed between the critical germ and the 
final plate size can be derived. Technically, an initially formed 
germ will be considered of size 2X2, randomly placed on the 
untransformed surface. Then it will ’’grow” to the final plate 
size, unless it encounters neighboring plates or the system bor¬ 
ders, case in which it only grows to the maximum size geomet¬ 
rically possible. 

Let us now analyze the consecutive transformations schemat¬ 
ically depicted in Fig.2 (a)-(c), which is just a more sophisti¬ 
cated version of Fig.l. On the first row, from left to right a 
direct transformation is simulated which consists in formation 
of plates with decreasing sizes as temperature is lowered, fi¬ 
nally resulting in the sizes distribution plotted in the bar chart. 
The distribution is balanced, in the sense that the areas of the 
product phase being in the form of different sizes squares is 
comparable (for instance there are 6 plates of side 12 with a to¬ 
tal area of 864, and 70 plates of side 3 with a total area of 630, 
etc.). 

The sizes distribution presented in the bar chart Fig.2a can for 
instance be compared with the chart distributions given in 1241 . 
to which a resemblance can be noticed for big size plates. How¬ 
ever, for small size plates ll24l gives a decrease not captured in 
our simulation most likely because: (i) 1241 addresses incipi¬ 
ent stages of the phase transition (usually few percents trans¬ 
formed, and mostly 30%) and (ii) the smallest plates size in our 
model was imposed by the discretization used, so the number 
of plates with the size ”3” should enclose also smaller plates 
from the physical situation. The size distribution of austenite 
nano-crystals is discussed in, e.g. ||25l. 

After the direct phase transition was simulated, we address 
the reverse transformation which is initiated with the smaller 
plates being the first to transform back and the system going 
basically through the same configurations in reverse order. In 
order to simulate the thermal arrest, this reverse transformation 
is incomplete. The sizes of the plates remained untransformed 
are shown in a rectangle. They define the configuration from 
which the next direct transformation will start (Fig.2b). One 
can notice an anomalously large number of big plates, con¬ 
sisting of the mentioned plates (remained untransformed) plus 
newly formed ones. An important consequence of this will be a 
depletion in the number of intermediate sizes plate, which can 
now fit with difficulty in the puzzle set by the large number of 
big plates. As such, most plates will grow to a lower size than 
the one maximum allowed at a given temperature. Repeating 
the ’’arrest” procedure enhances the effect (known as ’’hammer 
effect’ ’Eia), and this situation is shown in Fig.2c. 

The physical mechanism for depletion of intermediate sizes 
is the one schematically depicted in Fig.l, originating in the in¬ 
creased role of geometrical constraints in the direct transforma¬ 
tions following a thermal arrest. We present in Fig.2d the tem¬ 
perature dependence of the relative number of available places 
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Figure 2: (a) Left to right: a direct transformation resulting in the plates distribution sizes plotted in the bar chart. Next an incomplete reverse transformation is 
performed leaving the plates with sizes bigger than 8 untransformed -indicated with a rectangle in the bar chart-, (b) A subsequent direct transformation leads to 
a different plates distribution with increased number of big plates (sizes 9 to 12) and a depletion of intermediate sizes (sizes 7 and 8) . (c) Repeating the arrest 
procedure a second time accentuates the sizes unbalance, (d) Relative number of places where a square of size 8 can fit, for the three situations (a)-(c), as the 
temperature is lowered in the direct transformation (see description in text). Upper arrow indicates the transformation sense. 


to fit a size 8 square for the transformations (a)-(c). This num¬ 
ber is calculated by dividing the number of available spots to 
fit a size 8 plate with the number of places for the size 3 plates 
(smallest in our model). At the ignition of the initial transfor¬ 
mation (upper black curve) the ratio is nearly 1 , as all sizes can 
equally fit everywhere, except for near the system borders. As 
the transformation proceeds, smallest plates (size 3) fit signifi¬ 
cantly easier than larger plates (size 8 ) and the ratio decreases. 
After the arrest procedure, size 8 plates fit even more difficulty, 
which will determine them to grow to a smaller size, similar to 
the scheme shown in Fig.l. For all temperatures, the number 
of available places for the size 8 plates in transformations (b) 
and (c) are much reduced as compared to transformation (a), 
leading to the depletion seen in the bar charts. It is interesting 
to notice that the (a), (b) and (c) lines in Fig.2d start at progres¬ 
sively lower MS temperature, suggesting that that the thermal 
arrest may slightly increase the temperature hysteresis. 

In this section, a number of assumptions have been made 
regarding the direct transformation, however regarding the re¬ 
verse transformation it was only assumed that the plates trans¬ 
form back in the reverse order of their sizes. A possible trans¬ 
formation law example will be given in the next section. 

4. Reading thermal memory by reverse transformation (toy 

model illustration) 

Performing reverse transformations is a key part in the ar¬ 
rest experiments. In this section we assume a phenomenologi¬ 
cal transformation rate to simulate DSC curves for the reverse 
transformation. 


The simplest way to introduce the property (iii) is to assume 
that the interface energy has an elastic component, which can 


actually be responsible for initiation of the reverse product- 
parent transition. Let us denote the elastic fraction of the border 
energy with ALcr' (necessarily, one must have cr' < cr). This 
negative elastic energy that can be restored to the system will 
gradually make the martensite plates become energetically un¬ 
favorable, at a certain temperature, that is size dependent, and 
can be easily calculated to be: 


^ ^ yj Y v“/ 

eL 

Note that, indeed, the above formula implies automatically 
that the smaller martensite plates become unstable first (at lower 
temperatures) and will therefore be the first to transform back 
into austenite. Moreover, the reverse transformation is initiated 
below To, and possibly completed also below Tq. (For a dis¬ 
cussion of this possibility that the M ^ A transformation is 
initiated at A 5 ' < Tq, due to the stored elastic stress, see for in¬ 
stance (261). Once the condition Eq.2 is achieved for a given 
size plates, they begin to transform back with a certain transfor¬ 
mation rate, which can be generically expressed: 

/m^a(T, T) = J'^Nl 0(T - Tl)KL, n (3) 

where is a constant, Nl is the number of existing marten¬ 
site squares of side length L and 0(x) = 1 for x > 0, and 0 oth- 
ewise. The most simple choice for /(L, T) may be to consider 
it constant, alternatively one may use exponential form such as 
f{L, T) = exp[-^p-], with AQ(L, T) = -eL^{T-T q)-Act'L, 
as we employed for the sake of a graphical illustration in Fig.3. 
We stress again that Eq.3 fulfils the property (iii)- from Section 
2 , however we do not claim any microscopic justification for it, 
nor for the particular form of /. An important difference be¬ 
tween Eq.3 and Eq.l is that the M ^ A reverse transformation 
can actually stop if there exists no thermodynamically unstable 
squares of a given side L, and the transformation of the bigger 
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Figure 3: Simulated calorimetric signal (in arbitrary units) of the reverse transformation after previous complete cycles (a), after one previous thermal arrest at the 
temperature ”7^” (b) and after two successive arrests at the same temperature (c). 


squares cannot start until the temperature is further raised and 
the condition in Eq.2 is fulfilled. 

The numerical results are reflected in Fig.3, panel a) corre¬ 
sponding to a reverse transformation following complete previ¬ 
ous transformations. Panel b) corresponds to a reverse transfor¬ 
mation following an incomplete one (which was arrested at the 
temperature ’T^”) and a complete direct transformation. Panel 
c) illustrates the ’’hammer” effect, with a more pronounced dip 
after the arrest sequence is repeated. The plotted curves show 
also a number of oscillations, not related to the thermal mem¬ 
ory - the small local dips corresponding to the plates of a cer¬ 
tain size being transformed almost completely while the plates 
immediately smaller did not start the transformation. This is 
a limitation of the model, as we considered only 10 different 
sizes for the plates, a better discretization being perhaps neces¬ 
sary for quantitatively accurate plots. 


5. Relevance and limitations of the model 

A number of experimental papers (see, e.g. |[l0l[TTl) present 
microscopic evidence that, when increasing the temperature, 
the martensite plates transform back to austenite in reversed or¬ 
der of their formation. More general, even if some alloys un¬ 
dergo morphologic changes during cooling, the reverse trans¬ 
formation takes places through the same stages as the direct 
one - but in reverse order IZTll . Also in the paper {TH , chosen 
as example, the author comments on the fact that during the di¬ 
rect transformation, larger martensite plates are formed first -at 
higher temperatures- and smaller plates subsequently -at lower 
temperatures-. These aspects justify our assumptions described 
in Section 2. However, they are not universal features valid for 
all martensitic transformations. Also, one can in principle ar¬ 
gue that the fact that smaller plates form at lower temperatures, 
when the phase transformation is in an advanced stage, may in 
principle be exclusively due to geometrical constrictions and/or 
accumulated stress in the sample, rather than an intrinsic ther¬ 
mal property, as assumed in our numerical simulations. 

Further limitations of our numerical approach reside in the 
simplified shape choice (squares) of the plates, so more in¬ 
volved numerical approaches (in future works) may assume 
more realistic lamellar geometries and finite thickness. 

Even more drastically are the simplifications regarding the 
nucleation rates. While, for instance, it seems justified to 
assume isothermal martensitic transformation oa, the solid- 
solid nucleation rates are more complex then the simplest 
Eq.l assumed, with possible important consequences. Very 


recent state-of-the-art simulations actually suggest that solid- 
solid phase transitions take place in a two-step process, with an 
intermediate liquid-like phase Ga. 

Nevertheless, our intention here is just to give plausibility 
arguments for a mechanism potentially relevant for the temper¬ 
ature memory effect. 

6. Conclusions 

We propose a simple phenomenological model for a solid 
state phase transition, inspired by some experimentally ob¬ 
served features of the martensitic transformation. Then it is 
shown than the model naturally exhibits thermal memory. 

First, the key assumptions are described, namely that the di¬ 
rect phase transition takes place by formation of finite plates, 
whose seizes are temperature dependent (with bigger plates 
formed first) and also -possibly- geometrically restricted by ex¬ 
isting surrounding plates puzzle. Secondly, numerical simula¬ 
tions are performed considering isothermal nucleation rates for 
the direct transformation and the assumption that the reverse 
transformation takes place in reverse order (i.e. smaller plates 
transform back first). 

Based on the described model, if a reverse transformation 
is incomplete (arrested), the larger plates will remain untrans¬ 
formed and a subsequent direct transformation will show an 
anomalously large number of big sizes plates, consisting of the 
untransformed ones plus the newly formed. As temperature is 
further decreased, intermediate size plates may encounter dif¬ 
ficulties to geometrically fit - which results in them growing 
to a smaller final size and a general depletion of the sample 
of intermediate sizes. The effect is accentuated if the arrest 
scheme is repeated (the hammer effect). The dependence of 
the plates sizes distribution on the thermal history is a memory 
effect. Finally, a numerical model is proposed also for the re¬ 
verse transformation, showing that the memory effect should be 
’’readable” by a calorimetric signal. 

In a sense our approach has similarities with ||9j[T0|, the size 
of the plates playing the role of stress density, and each new 
direct transformation has a ’’stress free” start - or in our model 
big plates are formed first. However, the supplementary key 
role of geometrical constrictions imposed by the surrounding 
plates is a distinct feature. 

We hope that our phenomenological model may motivate fur¬ 
ther theoretical efforts to better understand the thermal memory 
effect, and also new experiments in particular seeking a corre¬ 
lation between incomplete thermal circles and the distribution 
of plates sizes. 
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